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A  Further  Generalized  Kettele  Algorithm 
With  Multiple  Constraints 

A.  Charnes 
D.  Zhang 

Abstract 

J.D.  Kettele  [1],  .using  dynamic  programming,  developed  a  simple  algorithm  (KA)  for 
the  optimal  redundancy  problem  in  reliability  and  life  testing  problems  with  a  single 

^  _  i 

constraint.  Proschan  ancTT.B;  Bray  \2]  gave  a  generalization  of  Kettele’s  dynamic 
programming  algorithm  to  include  multiple  constraints.  To  solve  a  much  broader  class  of 
optimization  problems  than  in  [1],  R.E.  Barlow  and  F?  Proschan  generalized  the  Kettele 
algorithm  (GKA)  to  apply  to  strictly  increasing  separable  function  problems  with  a  single 
constraint  [3]1 

In  this  paper,  we  consider  a  still  more  general  optimization  model  and  develop  a  Further 

..  -t 

Generalized  Kettele  Algorithm  (FGKA)  to  apply  to  multiple  constraints,  etc.  As  an  example,  an 
integer  Lexicographic  programming  model  will  be  solved,  (coronary  T,  section  2).^  Further¬ 
more,  another  form  of  the  more  general  optimization  model  is  pointed  out  in  section  4  of  the 
paper.  f  / 

\ 
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Kettele  [1]  presents  an  algorithm  for  allocating  redundancy  so  as  to  maximize  system 

reliability  without  exceeding  a  specified  cost.  Specifically,  a  system  consisting  of  k  "stages"  is 

considered.  The  system  functions  if  and  only  if  each  stage  functions.  Stage  i  consists  of  nj  (to  be 

determined)  units  of  type  i  in  parallel,  so  that  stage  functions  if  and  only  if  at  least  one  of  nj 

units  of  type  i  functions,  i  =  1,  2 . k.  Suppose  unit  i  has  a  "cost"  q,  i  =  1,...,k.  A 

unit  of  type  i  has  probability  p,  of  functioning,  independently  of  the  functioning  or 

nonfunctioning  of  the  other  units  of  the  system.  Thus  system  reliability  is  given  by 
p  (n)  =  n 1  =  ^  (!1  -(1  -Pi)n'j-  Then  the  problem  considered  in  [1]  is 

max  III1  'I1  ” P <)n ')• 

i-1 

s-«.|ciniSc.  (1) 

nj  >  0,  integer, 
where  c  is  a  limit  "cost”. 


Proschan  and  Bray  [2]  generalize  KA  to  solve  the  more  general  problem  of  maximizing 
system  reliability  without  exceeding  any  of  several  linear  constraints,  i.e.,  specifically, 

max  nl1  'I1  -Pi)"')- 

i-i 

k 

s.t.  ^CijniSCj,  (j- 1.2 . t) 

i-i 

nj  >  0,  integer, 

where  cy  is  the  "cost”  of  the  unit  i  of  the  jth  type,  i  =  1 . k;  j  -  1 ....  r,  and  Cj  is  the 

limit  "cost"  of  the  jth  type.  As  an  example,  the  first  type  of  cost  might  be  money,  the  second 
weight,  the  third  volume,  the  fourth  population  rate. 

Barlow  and  Proschan  [3]  then  give  a  Generalized  Kettele  Algorithm  to  solve  an 
optimization  model  more  general  than  in  [1]  as  follows.  Suppose  x-|, ....  x|<  are  k  variables 


called  "decision  variables";  x;  e  Sj  =  '<x\  '<  ...p,  i  =  1 . k.  Letxi  =  (x, . x,),  i  =  1, 


.  .  .  ,  k.  Assume  fj,  .  .  .  ,  f|<  are  strictly  increasing  functions,  with  yi(xi)  =  fi  (xi), 


y2(x^  =  f2  (yi ,  x2) ,  y3(x3)  =  f3(y2,  x3),  . . . ,  yk  (xk)  =  fk(yk.,,  xk).  Similarly,  assume  gi . Qk  are 

strictly  increasing  functions,  with  z^x,)  =  g^x,).  z2(x^  =  g2(z,,  x2), 
z  3(^3)  =  g3(z2,  x^, . . . ,  zk(xk)  =  gk{zk.1  ,xk).  It  is  suggestive  to  call  x,  an  "allocation"  of  order  i, 


yi  a  "payoff  of  order  i,  and  z;  a  "cost"  of  order  i  (i  =  1 . k)  even  though  the  model  is 


more  general  than  these  terms  would  indicate. 


The  Barlow  and  Proschan  (B-P)  model  is 
max  yk(xk), 


s.t.  zk(xk)<c, 


Xk=  (x, . xk)e  Six . xSk, 


where  c  is  a  limit  "cost".  In  this  model  there  is  only  a  single  constraint.  The  cost  Z|<  may  not 


be  a  sum  of  linear  functions  Cjxj ,  j  =  1,  .  .  .  ,  k.  It  is  strictly  increasing  in  the  number  of 


spares  of  each  type. 


B-P  gave  some  examples  (see  examples  and  exercises,  [3])  and  pointed  out  that  a  great 


many  other  models  requiring  optimization  subject  to  a  constraint  arising  in  reliability,  and 


more  generally  in  operations  research,  are  special  cases  of  the  general  optimization  model. 


Such  problems  may  be  solved  by  the  Generalized  Kettele  Algorithm  [3]. 


In  this  paper,  we  propose  a  more  general  optimization  model  and  develop  an  associated 


Further  Generalized  Kettele  Algorithm  (FGKA).  Thereby,  the  optimal  problems  which  can  be 


solved  by  FGKA  will  not  confine  to  a  simple  numerical  objective  function  or  a  single  constraint. 


Example  t:  An  Integer  Lexicographic  Problem 
max  (yk’{xk),  ylffxk), . yi'^Xk}), 


(i)  i- 

7-k  |xk  tSCi,  i  =  1,2 . m, 


M 


where  q  is  a  constant,  i  =  1,  . .  .  ,  m;  Xj  =  (xi . xj),  and  xj  is  a  real  variable, 

i,  j  =  1,  .  .  .  ,  k;  f VI ....  fk1  al1  are  strictly  increasing  functions  with  real  values,  and 

y i°(x -1}  =  f <J>(x  1) ,  yV' (x2}=f2J) (yV*. x2) . ykl)(xk)=fk1(yk-i.xk),  j=i . I;  gVl . . . , gL°  all  are 

(i),  ,  (j),  ,  (i)  (-  \  (i)  (  (i)  \ 

strictly  increasing  functions  with  real  values,  and  (Xi)  =  gi  (x -i) ,  z2  tx 2?  =  g 2  \z-\ 

ZkJ)  (xk)  =  gV)(zV-1 .  Xk),  z!-iUg<i1)(z[.i)1,x[1!),i>2,  j=1 . m;  (y^ . y['1  belongs  to  l-dimensional 

real  vector  space  R1  with  the  Lexicographic  ordering. 

The  Lexicographic  problem  is  a  special  case  of  our  still  more  general  optimiztion  model. 
And  as  we  know,  an  integer  bottleneck  problem  and  an  integer  time-cost  problem  are  in  some 
sense  special  cases  of  the  integer  Lexicographic  problems. 


Now,  let  us  give  the  still  more  general  optimization  model.  Suppose  Sj  is  a  countable 

{(1)  (2)  \ 

Xj  <  Xj  <  —  j,  i=l,2 . k.  Y,  Z  both  are  ordered  sets.  In  Z  there 

exists  a  maximal  element  z0  (or  °°).  Let  xj  represent  a  variable  called  a  "decision"  variable, 

xj  e  Sj,  and  xi  =  (x , , . . . Xj).  Let  fi  ,  f2 . fk!  Si .  92.  •  •  • .  9k  a11  be  strictly  increasing 

functions,  with  yi(xi)  -  fi(xi),  y2{x2)  =  f  2  (yi .  x2) . yk  (xk)  =  f  k  (y  k-i .  Xk}  andyi  e  V;  zi(xi) 

=  gi(xi),  z2(x4=92 (zi  , x2) . zJx^gklzk-LXk)  and  Zj  e  Z,  i  =  1,  2,  .  .  .,  k;  for  i  >  1 , 

gi(z0,  xj)  =  z0,  zj.1=  gj(zj_i ,  xi(1)). 

Call  x,  an  allocation  of  order  i,  yj  a  payoff  of  order  i,  and  Zj  a  cost  of  order  i, 
i  =  1,  2 _ k. 

The  still  more  general  optimization  model  is 
Max  yk (xk) , 

s.t. 

z.Msc,  (4) 

where  c  e  Z\{z0},  Xk^x,  ,  x2 . xk),  Xj  e  Sj  ,  j  =  1,  2 . k. 

Assume  that  (4)  and  Zj>  c,  j  =  1 , . .  .  ,  k  all  have  solutions. 


Allocation  xi  dominates  allocation  Xj’  if 
yi(xi)>yixi  )and  Zj(xi)< zixi ) < z0, or 
yi(xi)=yi(xi  )and  Zi(xi}< zj(xj  )<z0,or 


We  write  xj>Xi  .  We  also  say  the  corresponding  payoff-cost  pair  (yj,  Zj)  dominates 
(ys.  Zf)  and  write  (yi,  Zi)>(yi,  zj  . 


is  called  an  undominated  allocation  of  order  i  if  there  exists  no  Xj’  such  that 


A  _ 

Xi  >Xi  •  We  also  say  the  corresponding  payoff-cost  pair  (yj ,  Zj)  is  undominated. 


♦a'in  ,"^1 


‘I1 

:S 


•a 
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a 

•a 
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L  A  complete  sequence  of  undominated  allocations  of  order  i  (ending  in  x,  ,  say)  is 


-.(1)  -.(s) 

a  sequence  of  undominated  allocations  Xi . Xi  each  of  order  i,  such  that 


(JM 

(i) 

ydxi  )<■■ 

•  .<yi(xi  j 

(ii) 

ZilXi  }<.. 

■  .<Zi(xi  ) 

(iii)  If  xi  is  undominated  and  yields  a  payoff-cost  pair  distinct  from  those  of 

-I1*  -(s)  /->  Ms)\  (-i  Ps>\ 

xi  , . . . ,Xi  ,  then  yi(xi)>yj(xi  /  and  Zi(xi)> ZjVxi  ). 

I\  (1)  <1)\  I  (s)  (s)\ 

We  call  the  corresponding  sequence  (Vi  ,  Z\  ) . (yi  ,  2,  )  a  complete  undominated  sequence  of 


payoff-cost  pairs  of  order  i. 

We  now  present  an  algorithm  to  solve  this  more  general  optimization  model  (4).  We 

shall  call  it  the  Further  Generalized  Kettele  Algorithm  (FGKA). 

(i)  /  (i»  (0  /  (i)\  (si)  .  (s,+  i) 

Step  1 :  Compute  Vi  ■f1(x1J,z1  =  g i  ( x i  1  for  i  =  1 , 2,  .  .  .  ,  si ,  where  z,  <c,  and  z,  >c. 

There  exists  such  si  since  zi(xi)  =  Z|<(xi,  X2<1K  ....  x^1))  and  Z|<  <  c,  zi  >  c  both  have 
solutions.  (yV1,  zV’) . (y<i*,).z<1,1>)  constitute  a  complete  sequence  of  undominated  payoff-cost 

pairs  of  order  1  not  violating  the  cost  constraint. 

Step  2:  Compute  y2ij  =  f2  (yV).  ^L11) ,  z2ij=g2  (zV1.  x(2Jl)  for  i  =  1,  2 . si ;  j  =  1,  2,  .  .  .  such  that 

Z2ij  <  c.  Enter  the  payoff-cost  pair  (y2ij»  Z2ij)  in  row  i,  column  j  of  a  table  of  payoff-cost 


pairs. 


/  (i)  (i)\ 

Step  3:  The  lowest  cost  undominated  pair  (y2  ,z2 1  in  the  table  is  clearly  (y211l  z2n)  if  z211  < 
c.  If  z211  >  c,  there  is  no  solution. 

To  determine  [yf,  zf),  find  (y2ij .  z2ij)  such  that  yz^y^with  z2ijminimum  among 
entries  satisfying  this  inequality,  if  several  entries  of  identical  cost  qualify,  choose  the  highest 

payoff  y2ij.  If  several  entries  of  identical  cost  and  identical  highest  payoff  qualify,  choose  one 

/  (2)  (2)\ 

at  random,  since  all  are  equally  cost-effective.  The  payoff-cost  entry  so  chosen  is  (y2  . z2  j. 

I  M  (0\  /  (r  +  1)  (r  +  1)\  . 

Continue  in  this  fashion,  i.e.  having  found  (y2.z2/,then  (y2  ,z2  j  is  the  pair 

i  w 

(y2ij,z2ijJ  such  that  y2ij>  y2  with  z2ij  minimum  among  pairs  satisfying  this  inequality.  In  case 
of  ties,  follow  the  rules  for  breaking  ties  described  above. 


_  /  (®2>  (Sg)  (S2+1) 

Stop  at  \y2  .z2  /.where  s2  is  determined  so  that  z2  s  c,  while  z2  >c.  Such  S2 
exists  if  there  is  a  solution  to  (4). 

_  /  (V  h)\  f  (S  2)  M 

The  sequence  ^y2  .  z2 ^y2  .  z2  /  obtained  this  way  constitutes  a  complete  sequence 

of  undominated  payoff-cost  pairs  of  order  2  not  exceeding  the  cost  limit  c. 

Step  4:  Proceeding  in  a  similar  fashion  using  the  payoff-cost  pairs  (y3ij.  z3i|),  where 
.  /  (')  /  0)  (i)\  ,  .  _  ^ 

y3ij  =>3\y2  .z3j,  z3jj=g3\z2  ,x3j  for  1  =  1,  2,  .  .  .  ,  S2;  j  =  1 ,  2,  .  .  .  ,  such  that  z3ij  <  c,  obtain 
a  complete  sequence  of  payoff-cost  pairs  not  violating  the  cost  constraint.  Continue  in  this 
fashion  until  at  the  kth  stage,  arrive  at  the  complete  sequence  of  undominated  payoff-cost  pairs 
lyk.Zkj . \yk  ,zk  /,  where  zk  sc,  while  zk  >c. 

jg  \ 

Step  5:  Finally,  yk  is  the  maximum  payoff  achievalbe  under  the  cost  constraint;  the 

(sj  _,(**)  /  (s  k)  (s  k)\ 

corresponding  cost  is  zk  .  The  allocation  x  yielding  the  payoff-cost  pair  ^yk  ,zk  /  is  the 
solution  to  (4). 

Theorem  1.  The  payoff-cost  pair  (yk  k),zk  k))  and  the  allocation  obtained  using  the 
FGKA  is  the  solution  to  the  more  general  optimal  model  (4).  That  is,  it  is  the  maximum  payoff 

(S  k) 

achievable  with  corresponding  cost  zk  sc,  the  cost  constraint. 

Proof.  The  argument  of  the  proof  is  similar  to  that  of  the  original  theorem  about  the  GKA  (p. 
221,  [3]).  In  the  original  proof  two  lemmas  are  used.  They  are  as  follows. 

Lemma  4.5.  fat  Every  payoff-cost  pair  of  order  2  obtainsd  by  GKA  is  undnminated  (b) 
Every  undominated  payoff-cost  pair  of  order  2  may  be  obtained  by  the  GKA. 

LsmmaL.4,6.  Let  <y2,  z2)  $>  (y2\  z2’),  y3  =  f3  (y2,  x3),  z3  =  g3  (z2,  x3),  y3'  =  f3(y2\  x3), 
and  z3'  =  g3  (z2'f  x3).  Then  (y3,  z3)  £  (y3‘,  z3'). 


The  two  lemmas  are  correct  for  FGKA.  The  proof  of  the  former  is  almost  word  for  word 
repeat  of  the  original  lemma  proof.  The  proof  of  the  latter  is  little  different  from  the  original. 
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By  the  hypothesis  of  Lemma  4.6, 


(i)  y2  >  y2‘  and  z2  <  z2'  <  z0,  or 


(ii)  y2  =  y2'  and  z2<  z2’  <  z0,  or 


(iii)  z2'  =  z0. 


In  cases  (i)  and  (ii),  we  know  (y3,  z3)  %  (y3’,  z3’).  In  case  (iii),  this  also  holds  since  g3  is 


monotone  and  z3'  =  g3(z2',  x3)  =  z0.  Hence,  Lemma  4.6  holds  for  FGKA. 


Using  the  two  new  lemmas  and  repeating  the  argument  of  the  proof  of  the  original 


theorem  4.7  word  by  word,  we  get  a  proof  of  Theorem  2. 


Q.E.D. 


Now  we  show  that  the  general  optimization  model  (B-P)  under  multiple  constraints  is  a 


special  case  of  the  still  more  general  optimization  model  (4). 


’.  The  problem 


ykfck) , 


(i)  » 

zk  (xkJSCi,  1  =  1,2 . m, 

Xk=(x,,x2 . xkJeS1xS2x. .  .xSK, 


where  Sj,  yj,  xj,  Xj,  Yj,  all  are  as  in  the  more  general  optimization  model;  Z(')  is  an 


ordered  set  without  assuming  a  maximal  element,  is  a  special  case  of  the  still  more  general 


optimization  model  (4). 


Emflf.  Let  Zo  =  ZC)  x  Z<2>  x  ...  x  Z<m).  Define  an  equivalence  relation  =  among  elements  of  Z0  as 
follows.  For  a;,  bj  s  Z('),  i  =  1,  2 . m, 


(ai ,  a2,  .  .  . ,  am)  =  (b-j ,  b2,  . . . ,  bm) 


9 


if  and  only  if  either 

(i)  ai  >  Cj  for  some  i  =  2,  3,  .  .  .  ,  m,  and  bj  >  cj  for  some  j  =  2,  3 . m,  or 

(ii)  aj  <  Cj  and  bj  <  q,  i  =  2,  3 . m,  and  a,  =  b^ 

Thus  we  get  a  collection  Z  of  classes  of  equivalent  elements  of  Z0.  Denote  the  class  including  (clP 
c2 . cm)  by  c,  (zj’lxkj.zf’jxk) . zHxkljby zk(xk).  Define 

(alF  a2, - am)  >  (b1f  b2 . bm) 

if  and  only  if  either 

(i)  (a1t  a2 . am)  =  z0  *  (b1p  b2 . bm),  where  z0  is  the  class  whose  one  of 

components  2,  3,  ....  m  is  greater  than  corresponding  Cj,  or 

(ii)  z0  *  (a,,  a2 . am)  or  (b1t  b2 bm),  where  a ,  >  bv 

Then  Z  is  an  ordered  set  with  a  maximal  element  z0. 

Here  we  do  not  distinguish  a  clas"  from  its  representatives. 

So,  (5)  is  equivalent  to  a  (4)  in  which  zi<  e  Z,  c  e  Z\{z0}.  That  is,  (5)  is  a  special 
case  of  (4). 

Q.E.D. 


Similarly,  we  can  prove  corollary  1.  In  fact,  it  can  be  obtained  as  a  conclusion  of 
Theorem  1  itself,  also. 

Corollary  1.  The  integer  Lexicographic  problem  (3)  is  a  special  case  of  the  more 
general  optimization  model  (4),  or  (5). 


>.  Fill  Rate  Model  (compare  it  to  [3],  p.  209) 


IniCi<C0, 

i-i 

4 

Y  riiWi<  w0l 

i-i  (6) 

rij  >  0,  integer,  i  =  1,  2,  3,  4, 

where  n  =  (ni,  r\2,  n3,  n4)  is  a  spares  allocation;  c;  is  the  price  of  part  i,  c0  is  a  spares 
budget  constraint;  w;  is  the  weight  of  part  i,  w0  is  the  limit  weight  of  a  spare  parts  kit,  r(  n) 
is  the  fill  rate, 

Rlnl-tx/f  , 

1-1  j^>  I  ‘  (i-1  j 

where  s  0.  Here  we  assume  that 

(i)  Demands  for  spares  of  type  i  at  a  maintenance  depot  is  governed  by  a 

Poisson  process  with  demand  rate  Xj. 

(ii)  Enough  repair  facilities  are  available  at  the  depot  so  that  repair  of  a  failed 

unit  is  initiated  as  soon  as  it  is  received. 

(iii)  The  mean  time  to  repair  a  failed  unit  of  type  i  is  Uj. 

(iv)  The  cost  of  purchasing  n;  units  of  type  i  is  n-,cj. 

(v)  There  are  k  part  types,  i.e.,  i-1,2 . k. 

In  (6)  there  are  two  constraints,  which  says,  that  the  budget  is  not  to  be  exceeded  nor  is  the 
allowed  total  weight  (or  volume). 


To  solve  the  example,  first  let  Si  =  {0,  1, 2,  .  . . },  Y  the  set  R1  of  real  numbers.  To 
define  Z,  first  define  an  equivalence  relation  =  in  the  2-dimensional  real  vector  space  R2  as 
follows. 

(Ui,  U2)  s  (v1t  v2) 
if  and  only  if  either 

(i)  u1  >  w0  and  v,  >  w0,  or 

(ii)  u1(  v,  <  w0  and  u2  =  v2 

Denote  (u1t  u2)  by  z0  if  u,  >  w0.  Then  let  Z  be  the  collection  of  equivalence  classes  of  elements 
in  R2.  Define  an  order  relation  >  in  Z: 


( u i ,  u2)  >  (v1f  v2) 
if  and  only  if  either 

(i)  U!  >  w0  and  v,  ^  w0,  or 


(ii)  Ui ,  v1  ^  w0  and  u2  >  v2. 

Denote  (w0l  cQ)  by  c. 

The  decision  variables  are  to  be  the  number  of  spares  of  the  k  different  types,  (n1(  n2, 
n3,  n4).  The  functions  of  fj  and  gj  of  the  more  general  model  are  given  by 

j-0  J! 


MR.niJ-R  +  X.X^i-  e'"1'5’ 


i-o 


fli(ni)  =  (w1  n1t  Ci  n-i), 

9i  ((Wi  c)),  Hi)  =  (w,  c)  +  (Wj  nj,  Ci  nj). 


So,  all  fj  are  strictly  increasing  functions.  And  all  gj  are  increasing  in  Z,  strictly  increasing 
in  Z\(z0}. 
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For  example,  consider  the  situation  of  a  budget  c0  =  1500  (in  dollars),  weight  limit  w0 
=  100  (in  kg),  and  demand,  repair,  weight  and  cost  data  as  shown  in  Table  1. 


Part 
Type,  i 

Demand  Rate 
(per  hour) 

Mean  Time 
to  Repair 
Failed  Part 
vi 

kg  weight 
of  part, 

Wi 

Dollar  Cost 
of  part, 

Cj 

computed 
value:  l)j 

1 

.01 

100 

30 

200 

1.0 

2 

.02 

150 

10 

100 

3.0 

3 

.03 

60 

40 

300 

1.8 

4 

.01 

200 

5 

250 

2.0 

Table  1 .  Input  Data  for  Spares  Allocation 


K  — * 

To  simplify  the  computation  throughout,  we  drop  the  denominator  to  £i  ^  in  R(n)  since 


it  is  constant  throughout  and  just  compute  the  fill  rate  numerator  where 


1111  [X.v)1  -mi 

Ri(n)  =  Il- -Le  i  =  1 ,  2 . 


i-o 


The  steps  taken  in  carrying  out  FGKA  are  the  following. 


(1)  Noting  that  nj  =  0,  1,  2,  ....  in  Step  1  we  get 


l^e-"  Voicin'!  e'1 


i-0 


i! 


i-0 


and 


g i ( n i )  =  (w,  n,  ,  n,)  =  (30n1f  200n1), 


and 


c  =  (w0,  c0)  =  (100,  1500) 


1  3 


we  get  the  complete  sequence  of  undominated  payoff-cost  pairs  of  order  1 , 

(y  1 (1),  z1d))=  (0,  (0,  0)),  (yi(2),  z^2))  =  (.00368,  (30,  200)), 

(yi<3),  2lO))  =  (.00736,  (60,  400)),  (y i  <4>,  z<\  <4>)  =  (.00920,  (90,  600)). 
r-itsie  3-|  =4 

(2)  Noting  that  n2  =  0,  1,  2 . in  Step  2  and  Step  3, 


.  V  O.  3.0  J 

f2(Ri,  n2)  =  R  +  X2  S  — 7T~  e  =R  +  .02^-^-e 

i-o  ‘ ' 


i-o 


J! 


and 


g2((w,c),  n2)  =  (w,  c)  +  (10  n2  ,  100  n2), 

we  get  Table  2. 

The  complete  sequence  of  undominated  payoff-cost  pairs  of  order  2  is  shown  by 
- •  in  Table  2.  Here  s2  =  9. 


• - • - •- 


(3)  Since  n3  =  0,  1 ,  2,  .  .  .  ,  in  Step  4 


f3(R.  n3) 


"x1  (tauJ  n’l1  (l  fl|^  -1.1 

=  R  +  Xj  X  Hr-e  =  R+03X  ^  e 

j_0  1- 


j-0 


and 


g3((w,  c),  n3)  =  (w,  c)  +  (40  n3  ,  300  n3), 

we  obtain  Table  3.  Here  s<3>  =  9,  i.e.,  the  complete  sequence  of  undominated  payoff-cost  pairs  of 
order  3  shown  in  Table  3  has  9  elements. 


(4)  Noting  that  n4  »  0,  1,  2,  ....  in  Step  4 


n4,1  (X-4  U4) 


f4(R.  n4»  -  R  +  X.4  =  R  +  -01X 


n*'1  to 


(2.0)  _  —2.0 


j-0 


i-o 


j! 


and 


■mensnnwnBWWiiviwnvMviwiiviivi  vmu" 
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The  complete  sequence  of  undominated  payoff-cost  pairs  of  order  4  is 


o 

o 

o 

o 

0, 

1, 

0, 

0; 

0, 

2, 

0, 

0: 

0, 

3, 

0,  0; 

o 

o 

o 

o, 

5, 

0, 

0; 

1. 

4, 

0, 

0: 

1. 

5, 

0,  0; 

2,  4,  0,  0; 

2, 

4, 

0, 

1; 

2, 

4, 

0, 

2; 

1, 

5, 

0,  3. 

(5)  Finally,  a  solution  of  example  3  is 
( n i ,  n2,  n3,  n4)  =  (1,5,  0,  3). 
with  Max  (n)  =  .02576. 


(i) 

n  „ 


NS 

FRN 

WC 

0 

0 

(0,  0) 

1 

.00010 

(10,100) 

2 

.00548 

(20,200) 

3 

.00996 

(30,300) 

4 

.01332 

(40,400) 

5 

.01533 

(50,500) 

6 

.01634 

(60,600) 

7 

.01677 

(70,700) 

8 

.01694 

(80,800) 

9 

.01699 

(90,900) 

10 

.01701 

(100,1000) 

0 

0 

(0, 0) 

0,  0 

0,1, 

°'2i 

0,3t 

0.4  f 

0.5, 

0.6 

.01634 

(60,600) 

0,7 

.01677 

(70,700) 

0,8 

.01694 

(80,800) 

0,9 

.01699 

(90,900) 

0,  10 

.01701 

(100,1000) 

0 

(0,0) 

.00010 

(10,100) 

.00548 

(20,200) 

.00996 

(30,300) 

.01332 

(40,400) 

.0/3 

JM.SOO) 

1 

.00368 

(30,200) 

1,0 

.00368 

(30,200) 

1 ,1 

.00468 

(40,300) 

1,2 

.00916 

(50,400) 

1.3 

.01304 

(60,500) 

1.V£ 

.01700 

(70,600) 

.ojpfi 

JftoJOO) 

1,6 

.02002 

(90,800) 

1,7 

.02045 

(100,900) 

2 

.00736 

(60,400) 

2,  0 

.00736 

(60,400) 

2,  1 

.00834 

(70,500) 

2,2 

.01284 

(80,600) 

2,3 

.01732 

(90,700) 

2,4/ 

.02068 

(100,800) 

3 

.00920 

(90,600) 

3,0 

.00920 

(90,600) 

3,1 

.01020 

(1 00,700) 

Table  2  (y2ij,  *2ij) 


Notation:  NS  -  number  of  spares,  FRN  -  fil!  rate  numerator,  WC  -  weight  and  cost. 


Corresponding  to  (1),  consider  an  optimal  design  problem 

k 

min  VC|ni 

st.  n(i  -(*•  -Pi)n)s  r. 


nj  >  0,  integer, 

where  R  is  a  limit  reliability,  q,  nj,  pi  are  as  in  (1). 

Generally,  let  us  consider  a  resolvable  problem 
min  yk(xi), 

s.t.  zk  (xi)  >  c,  (  8  ) 

where  xk,  yk,  zk  and  c  are  as  in  (4).  The  only  difference  is  that  in  Z  there  exists  a  minimal 

element  z0,  but  a  maximal  element  Zj<  c  has  a  solution,  j  =  1 . k.  And 

Sj-{...<Xj(2)<xj(1)},  1,2 . k.  (9) 

For  (7),  select  positive  integers  mi . m*  large  enough.  Let  Sj  «  {  0  <  .  .  .  <  mj  - 1 

<  mj},  y-i(xi)  =  f i ( x t )  =  c^l  y2(x4  =  f2(yi.  x2)  =  yi  +  C2X2 . yk{xi)  =  fk(yk.i,  xj  = 

yk-i  +  ck  xk;  z^xO  =  (1  -  (1  -  pO*1),  z2  =  zt(1  -  (1-  p2)Xz)  —■  . 

(1  -(i-p2) 

1 

z k  =  zk.i(1  -  (1  -  pk)*k).  ..  mk.  Then  (7)  is  an  example  of  (8). 

(1  -(l-Pk) 

We  point  out  that  (8)  is  another  form  of  the  more  general  optimization  model  (4).  Let 
Sj*  =  {xj(1 )  zl  xj(2)  zi  . . .  }  be  an  ordered  set  with  ordering  A,  i  =  1,  .  .  . ,  k.  That  is  xj(u)  zi 
xj(v)  if  and  only  if  Xj(u)  >  x,(v).  Define  a  new  ordering  zl  in  Y  such  that  azib  if  and  only  if  a  > 
b.  Denote  the  Y  with  the  new  relation  zl  by  Y*.  Similarly,  we  get  Z*.  Thus  z0  is  a  maximum 
in  Z*  about  the  relation  zl.  Now,  the  problem  (8)  gets  a  new  form  as  follows, 
max  yk  (xk), 

s.t.  zk(xk)zi_c,  (10) 

where  the  sign  max  is  about  the  new  relation  zl  in  Y*  and  the  equality  sign  is  the  same  with  the 
original. 
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Theorem  3.  The  problem  (8)  is  just  another  form  of  the  still  more  general  optimization  model 
(4). 

Hence,  we  can  use  the  FGKA  solving  (8)  or  (10). 

Furthermore,  let  us  look  at  another  optimal  design  problem 

k 

min  Xc'n' 

i.i 


s.t. 


n(i  -{i-Pi)n')^R 


]T  w;n;  <w, 


(11) 


i.i 


where  wj  is  the  weight  of  a  part  of  type  i,  w  is  the  total  weight  limit.  Other  signs  are  as  in 

(7). 

Generally,  we  would  consider  a  problem 
min  yk(xk), 

S.t.  zk(xk)>C  (12) 

zk'  (xk)  £  c\ 

where  xk,  yk,  zk,  c  are  as  in  (8),  zk'  e  Z',  an  ordered  set  with  a  maximal  element  z0',  and 
c‘  e  Z'  -  {z0  } . 

Theorem  4.  The  problem 


min  yk  xk 


s.t.  zk 


(Of- 


z.(i) 

z  k 


(xk)^Ci  , 
(xk)^C'i 


i  =  1, 
i  =  1, 


m, 

n, 


(13) 


Xk=(x, . xk)e  Six . . .  xSk, 


where  Sj,  Y,  xj,  xif  yj,  all  are  as  in  theorem  2,  z*Jl>  is  as  z  j0  in  theorem  2.  The 
corresponding  gj’  takes  values  in  the  corresponding  set  Z’;  Zj(')  (xj)  =  gj(')(Zj.i  (•),  Xj)  is 
increasing  and  takes  values  in  an  ordered  set  Z  with  minimum  z0,  is  a  special  case  of  (12). 


Proof.  Let  Zq  =  Z^)x  . . .  xZ(m).  Define  an  equivalence  relation  =  among  elements  of  Z0  as 
follows.  For  aj,  bj,  e  Z('),  i  =  i . m,  (ai ,  .  .  .  ,  am)  =  (bi ,  .  .  .  bm)  if  and  only  if  either 

(i)  aj  <  Cj  for  some  i  =  2 . m  and  bj  <  Cj  for  some  j  =  2,  .  .  .  ,  m,  or 

(ii)  aj  >  Cj  and  bj  >  c;,  i  =  2 . m,  and  ai  =  bi . 

Thus  we  get  a  collection  Z  of  classes  of  equivalent  elements  of  Z0.  Denote  the  class  including 

(ci .....  Cm)  by  c,  (zkO  )(xk), .  . .  zk(m){xkj)  by  zk  (xk).  Define  (ai ,  . .  .  am)  <  (bi . bm)  if 

and  only  if  either 

(i)  (ai . am)  =  z0  *  (bi . bm),  where  z0  is  the  class  whose  one  of  components 

2,  3 . m  is  less  than  corresponding  q,  or 

(ii)  z0  *  (ai . am)  or  (bi,  .  .  .  ,  bm),  where  ai  <  bi. 

Then  Z  is  an  ordered  set  with  a  minimum  z0.  Similarly,  let  Z0’  =  Z'COx  .  .  .  xZ'(n).  Define  an 
equivalence  relation  s  among  elements  of  Z0'  as  follows.  For  aj,  bj  e  Z’('),  i  =  1, . . .  ,  n, 

(ai ....  an)  =  (bi ...  .  bn)  if  and  only  if  either 

(i)  a,  >  Cj’  for  some  i  ■  2, . .  . ,  n  and  bj  >  Cj'  for  some  \  -  2 . n,  or 

(ii)  a;  5  Cj',  bj  <  q',  i  =  2,  .  .  .  ,  n,  and  ai  =  bi . 

Hence  we  get  a  collection  Z '  of  classes  of  equivalent  elements  of  Z0‘.  Define  the  class  including 

(ci’ . cn)  by  c'  and  (zk' C )(xk), . . .  zk'(n)(xk))  by  zk'(xk). 

Define  (ai . am)  >  (bi, .  .  . ,  bn)  if  and  only  if  either 

(i)  (ai , ....  an)  =  z0'  *  (bi , ... ,  bn),  where  z0'  is  the  class  whose  one  of  components 

2,  ....  n  is  greather  than  corresponding  q',  or 

(ii)  zQ'  *  (ai . an)  or  (bi , .  . . ,  bn),  where  ai  >  bi . 

Then  Z'  is  an  ordered  set  with  a  maximal  element  z0'. 

Therefore  (13)  is  equivalent  to  a  (12)  in  which  Z|<  e  Z,  c  e  Z\(z0},  z*'  e  Z', 
c'  €  Z\{z0'}.  That  is,  (13)  is  a  special  case  of  (12). 
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